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. Abstract. Formulae suitable for computing LID acceleration and correspond- 

ing diffusional segregation of isotopes in CP star model atmospheres are given. 
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Light-induced drift (LID) as a phenomenon responsible for separation of isotopes 
C"| | in the atmospheres of CP stars has been proposed by Atutov and Shalagin 

£Li (1988). Thereafter we have studied evolutionary abundance changes of Hg and 

q its isotopes (see Sapar et ai, C_Sapar, these proceedings). Here we give formulae 

i ■ describing diffusion of isotopes in the form best suited for numerical model 

r/3 | computations. 

Diffusive transfer of radiation holds in the deeper layers of stellar atmo- 
spheres, where monochromatic radiative flux can be described by formula 

;> F v = FKr , where = / dv. 

\q • HwdT K R J K v dT 

lO ■ Here F is the total radiation flux and Kr is the Rosseland opacity integral. 

In the opacity coefficient n v the overlapping spectral lines of the trace element 
studied play an important role, and thus k v = c v + J2j a jWj(u„), where c„ 
is the continuous opacity coefficient and o~j is the transition cross-section per 
gram. Summation is made over spectral lines j with line profile functions Wj, 
usually being the Voigt functions with argument given by u v = (y — Vj)j Avt, 
where Apt is the thermal Doppler width of the spectral line. As we have shown, 
the effective acceleration of LID due to spectral line j can be expressed quite 
, similarly to the usual expression of radiative acceleration 

a. j = — / — ^ — -F v dv, % = qeo-j, (1) 
c Jo ou u 

where q = Mvrc/2hv and the efficiency of LID is e — (C u — Ci)/(A U + C u ). 
Here C u and Ci are the collision rates of particles in upper and lower states, 
respectively, and A u is the probability of spontaneous transitions. 

Expression (1) in the region of diffusive transfer of radiation takes the form 

J c J ou v n u dT 

Both isotopic and hyperfine splitting of spectral lines of all ions should be taken 
into account to calculate accelerations producing segregation of isotopes. 
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The equation of continuity for isotope i in the plane-parallel stellar atmo- 
sphere has the form dpi/dt + d(piVi)/dr — 0. The model atmosphere data cor- 
respond to standard points, being equidistant on the logarithmic scale of mean 
optical depth. These points are enumerated as layers n growing downwards. 
Treating n as a continuous parameter, we change variables in equation of con- 
tinuity. Since ^ = %lfdn an( ^ Pin ~ ~5n> wnere I 1 i s total column density, we 
get for radial gradient 

d d p p 

= —7—, where 7 = 



dr dn' dp/dn pd\np/dn 

Denoting a ratio of current concentration to its initial value as Cj, we can write 
Pi = p\C{ and the equation of continuity reduces to 

dlnCj = 7 djpjVi) 
dt pi dn 

Logarithms are used to avoid possible negative values of C, in time integration. 
Diffusion velocity V, in the presence of stellar wind can be found from equation 

T , , x, .dpi dpi dpi 
p i V i =p i (a i -g)t-A—, ^r = -7^, (4) 

where at is the sum of radiative and LID accelerations, t is the mean free flight 
time of the particles , m is the mean mass of buffer particles and A = kTt/m 
is the diffusion coefficient of trace particles. Thus from equation (4) we find 

Vi _ mjaj - g) ^ dhi(p"Ci) Vj_ _ ma, ^ dlnC'j 

7A kTj dn 7A kTj dn 

Here the last expression has been obtained from the first taking into account 
that for model stellar atmospheres approximately holds mg/kT-f = dlnp^/dn. 
The velocity Vi is to be used in the equation of continuity (3), reducing it to 
a generalized Fokker-Planck equation. The equation (5) can be used also for a 
crude prediction of concentrations d . For the final equilibrium state in the case 
of no stellar wind Vi — and in the case of constant mass loss rate piV, = rhi. 

As we see, for evolutionary computations we need to find derivatives dhip/dn, 
dp/dn and d^/dn, which correspond to buffer gases and several derivatives for 
each isotope of the trace (impurity) particles. The derivatives have been found 
using the 4th order Lagrange interpolation formulae for equidistant nodes. 
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